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Abstract 

We generalize the f(R) type gravity models by assuming that the gravitational Lagrangian 
is given by an arbitrary function of the Ricci scalar R and of the matter Lagrangian L m . We 
obtain the gravitational field equations in the metric formalism, as well as the equations of motion 
for test particles, which follow from the covariant divergence of the energy-momentum tensor. 
Lhe equations of motion for test particles can also be derived from a variational principle in the 

co : 

particular case in which the Lagrangian density of the matter is an arbitrary function of the energy- 
density of the matter only. Generally, the motion is non-geodesic, and takes place in the presence 
00 . . . ^ ^ 

of an extra force orthogonal to the four- velocity. The Newtonian limit of the equation of motion 



cr 



o 

is also considered, and a procedure for obtaining the energy-momentum tensor of the matter is 
presented. The gravitational field equations and the equations of motion for a particular model in 

X, 

which the action of the gravitational field has an exponential dependence on the standard general 
relativistic Hilbert-Einstein Lagrange density are also derived. 
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I. INTRODUCTION 

A very promising way to explain the recent observational data of the late-time 

acceleration of the Universe and of dark matter is to assume that at large scales the Einstein 
gravity model of general relativity breaks down, and a more general action describes the 
gravitational field. Theoretical models in which the standard Einstein-Hilbert action is 
replaced by an arbitrary function of the Ricci scalar R, first proposed in [3j, have been 
extensively investigated lately. Thus, the late-time cosmic acceleration can be explained 
by f{R) gravity |4|, and the conditions of viable cosmological models have been derived 
in |5|. In the context of the Solar System regime, severe weak field constraints seem to 
rule out most of the models proposed so far ja Q], although viable models do exist js-ll]. 
f(R ) models that pass local tests and unify inflation with dark energy were considered in 



121 ] . The possibility that the galactic dynamic of massive test particles can be understood 



without the need for dark matter was also considered in the framework of f{R) gravity 



models 



13 



17|. For a review of f(R) generalized gravity models sec |18]. 



A generalization of the f(R) gravity theories was proposed in |l9] by including in the 
theory an explicit coupling of an arbitrary function of the Ricci scalar R with the matter 
Lagrangian density L m . As a result of the coupling the motion of the massive particles is 
non-geodesic, and an extra force, orthogonal to the four-velocity, arises. The connections 
with MOND and the Pioneer anomaly were also explored. This model was extended to 
the case of the arbitrary couplings in both geometry and matter in [20]. The astrophysical 
and cosmological implications of the non-minimal coupling matter-geometry coupling were 



extensively investigated in 



21 



22] 



It is the purpose of the present paper to consider a maximal extension of the Hilbert- 
Einstein action, by assuming that the gravitational field can be described by field equations 
that can be derived from a variational principle with the Lagrangian density of the field 
given by an arbitrary function of the Ricci scalar and of the matter Lagrangian density. We 
derive the field equations for this model, and obtain the equations of motion of the test 
particles. Generally, the covariant divergence of the energy-momentum tensor is non-zero, 
and the motion of test particles is non-geodesic, and takes place in the presence of an extra 
force, orthogonal to the four-velocity. The non-geodesic motion, due to the non-minimal 
couplings present in the model, imply the violation of the equivalence principle, which is 
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highly constrained by solar system experimental tests 23|, |24|. However, it has been recently 



reported, from data of the Abell Cluster A586, that interaction of dark matter and dark 



energy does imply the violation of the equivalence principle [25|. Thus, it is possible to test 
these models with non-minimal couplings in the context of the violation of the equivalence 
principle. It is also important to emphasize that the violation of the equivalence principle 
is also found as a low-energy feature of some compactified version of higher- dimensional 
theories. 

The Newtonian limit of the equation of motion is also investigated, and a procedure 
for obtaining the energy-momentum tensor of the matter is presented. The gravitational 
field equations and the equations of motion for a particular model in which the action of 
the gravitational field has an exponential dependence on the standard general relativistic 
Hilbert-Einstein Lagrange density, f(R,L m ) = Aexp (R/2A + L m /A), are also derived. 
Generally, we assume that our model is valid for all types of matter, including scalar fields, 
vector fields, or even Yang-Mills fields. Non-minimal models of f{R) gravity, in which the 
geometric part of the action is multiplied by the electromagnetic field or the Yang-Mills field 



Lagrangian, were considered in 26]. 



The present paper is organized as follows. The field equations of the / (R, L m ) gravity 
model are derived in Section [TXJ The equation of motion of test particles, and the Newtonian 
limit of the equation of motion, are considered in Section [TTT1 The construction of the energy- 
momentum tensor for perfect thermodynamics fluids is presented in Section HVl A particular 
model of / (R, L m ) gravity, in which the Lagrangian density of the gravitational field has 
an exponential dependence on the standard Hilbert-Einstein Lagrangian density of general 
relativity is discussed in Section [V] We conclude and discuss our results in Section [VTl 

II. FIELD EQUATIONS OF f(R,L m ) GRAVITY 

We assume that the action for the modified theories of gravity takes the following form 



S 



J f(R,L m )^d 4 x , (1) 



where / (R, L m ) is an arbitrary function of the Ricci scalar R, and of the Lagrangian density 



corresponding to matter, L m . We define the energy- momentum tensor of the matter as 



27] 



By assuming that the Lagrangian density L m of the matter depends only on the metric 
tensor components g^, and not on its derivatives, we obtain 



(3) 



By varying the action S of the gravitational field with respect to the metric tensor com- 
ponents g^ v we obtain first 

5S =-. i f l{ (jR.L m ),)n + f l .. {n.L lll )^8g^ - 



:9nuf (R, L m ) Sg 



-gd 4 x, 



(4) 



where we have denoted f R (R, L m ) = df (R, L m ) jdR and f Lm (R, L m ) = df (R, L m ) /dL m , 
respectively. For the variation of the Ricci scalar we obtain 

5R = 5 {g^R pv ) = R pv 5g^ + <T (v^I^ - VJT* X ) , (5) 

where Va is the covariant derivative with respect to the symmetric connection T associated 
to the metric g. Taking into account that the variation of the Christoffel symbols can be 
written as 



6T 



A 



g a (VpSgva + VJg a/1 - V a 5g 
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we obtain for the variation of the Ricci scalar the expression 

5R = R^bgv + ^V M V^<T - V^vSgr. 
Therefore for the variation of the action of the gravitational field we obtain 
5S = J [f R (R,L m )R^5g^ + 



(6) 



(7) 



f R (R,L m )g^V^5g 



f R (R,L m )V^Jg^ + f Lm (R, L r 



SL r . 



Sg 



pi/ 



flV 



:9iiuf (R, L r 



-gd 4 x . 



(8) 



After partially integrating the second and third terms in Eq. (jSJ), and with the use of the 
definition of the energy-momentum tensor, given by Eq. ([3]), we obtain the field equations 
of the / (R, L m ) gravity model as 

f R {R, L m ) R p „ + ^V M V - V^V.) f R {R, L m ) - 
o [/ (R, L m ) - fL m (R, L m ) L m ] g^y = 



/lm (R> Lm) Tyiw 



(9) 



If f(R,L m ) = R/2 + L m (the Hilbert-Einstein Lagrangian), from Eq. (Q we recover 
the standard equations of general relativity, R^ u — {l/2)g llv R = T^ u . For f(R,L m ) = 
fi(R) + f2(R)G (L m ), where fi, f 2 and G are arbitrary functions of the Ricci scalar and of 
the matter Lagrangian density, respectively, we reobtain the field equations of the modified 



gravity with arbitrary geometry- matter coupling, considered in 20] . 

The contraction of Eq. ([9]) provides the following relation between the Ricci scalar R, the 
matter Lagrange density L m , and the trace T = of the energy- momentum tensor, 



f R (R, L m ) R + 3V^f R (R, L m ) - 2 [/ {R, L m ) - 

f Lm (R, L m ) L m ] = l -f Lm (R, L rn ) T. (10) 

By eliminating the term V^V^/r (R, L m ) between Eqs. (jSJ) and (JTOl . we obtain another 
form of the gravitational field equations as 

fa {R, L m ) (r^ — -Rg^J) + 

g [/ {R, L m ) - f Lm (R, L m ) L m ] g^y = 

-fh m (R, L m ) (t^ - ^ T 9fii^j + 

V,S?MR,L m ). (11) 
By taking the covariant divergence of Eq. ([9]), with the use of the mathematical identity 



28|, 



/r (R, L m ) R pv — -f (R, L m ) g^ v + 



(^V M V - V M V„) f R {R, L m )) = , (12) 



we obtain for the divergence of the energy-momentum tensor T^ v the equation 

V% = In [f Lm (R, L m )} {L m9tiU - T, u } = 

2W]n[f Lm (R,L m )}-^. (13) 

The requirement of the conservation of the energy-momentum tensor of matter, V^T^ = 
0, gives an effective functional relation between the matter Lagrangian density and the 
function f Lm (R, L m ), 

Wln[f Lm (R,L m ))-^ = 0. (14) 



Thus, once the matter Lagrangian density is known, by an appropriate choice of the 
function / (R, L m ) one can construct, at least in principle, conservative models with arbitrary 
matter-geometry dependence. 



III. THE EQUATION OF MOTION OF TEST PARTICLES AND THE NEWTO- 
NIAN LIMIT IN f(R,L m ) GRAVITY 



As a specific example of / (R, L m ) gravity model we consider the case in which the matter, 
assumed to be a perfect thermodynamic fluid, obeys a barotropic equation of state, with the 
thermodynamic pressure p being a function of the rest mass density of the matter (for short: 
matter density) p only, so that p = p(p). In this case, the matter Lagrangian density, which 
in the general case could be a function of both density and pressure, L m = L m (p,p), or of 
only one of the thermodynamic parameters, becomes an arbitrary function of the density of 
the matter p only, so that L m = L m (p). Then the energy-momentum tensor of the matter 
is given by [20] 



T v> = p d J^ u ^ + ( Lr 



' v- p ^) r - (15) 

where the four-velocity = dx^/ds satisfies the condition g^u^Uy = 1. To obtain Eq. (fl5l) 
we have imposed the condition of the conservation of the matter current, V„ {pu u ) = 0, and 
we have used the relation 

(16) 



whose proof is given in the Appendix of 22J. With the use of the identity 

d 2 x fl 



ds 2 



(17) 



from Eqs. ( TL3l) and ( IT5l) we obtain the equation of motion of a test fluid in / (R, L m ) gravity 
as 

+ rVV = r, (is) 

where 



- V„ In 



fo m (-ft) L r 



dL m (p) 
dp 



-g 
o. 



/a/ \ 



(19) 



The extra-force is perpendicular to the four-velocity, f^u^ 

The equation of motion Eq. (fI5|) can be obtained from the variational principle 



5S P = S J L p ds = S J \fQ\J g l _ iu u^u u ds = 



(20) 



where S p and L p = \fQ^/g^ v u^u v are the action and the Lagrangian density for test particles, 
respectively, and 

Q = f Lm {R,L m ) d -^. (21) 
To prove this result we start with the Lagrange equations corresponding to the action (1201) . 



d ( dL p \ dL r 



ds \du x J dx } 

Since 



0. (22) 



du x 
and 

dL 



° L " 'Qu x (23) 



dx x 2 V ^ M "' A 2 Q 

a straightforward calculation gives the equations of motion of the particle as 



f/s , + K x u u u x + (u»u v - gT) V, In yjQ = 0. (25) 

By simple identification with the equation of motion given by Eq. ( ITS]) , we obtain the 
explicit form of \/Q, as given by Eq. (j2T|) . When \/Q — > 1 we reobtain the standard general 
relativistic equation for geodesic motion. 

The variational principle (j2U|) can be used to study the Newtonian limit of the equations 
of motion of the test particles. In the limit of the weak gravitational fields, 

ds RJ ^l + 2(/)- iPdt w (l + (j) - v 2 /2) dt , (26) 

where <ft is the Newtonian potential and v is the usual tridimensional velocity of the fluid. By 
assuming that in the Newtonian limit of weak gravitational fields the function fi m (R, L m ) 
can be represented as 

f Lm (R,L m )^l + U(R,L m ) , (27) 
where U (R, L m ) 1, the function y/Q can be obtained as 

dp dp 
Taking into account the first order of approximation, the equations of motion of the fluid 
can be derived from the variational principle 
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dL m {p) dL m (p) v 2 

^p- + f/(/? ' Lm) ^p- +0 -y 



dt = 0, (29) 
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and are given by 



dp 



(30) 



where a is the total acceleration of the system, 
acceleration, and 



a E = -VU E = -V 



U(R,L m ) 



-V0 is the Newtonian gravitational 
dL m (p) 



dp 



(31) 



is a supplementary acceleration induced due to the modification of the action of the grav- 
itational field. As for the term an = —V[dL m (p) /dp], it has to be identified with the 
hydrodynamic acceleration term in the perfect fluid Euler equation. 



IV. THE MATTER ENERGY-MOMENTUM TENSOR IN / (R, L m ) GRAVITY 

Since % can be interpreted as a hydrodynamic acceleration term, in the Newtonian limit 
it can be represented in the form 

a, = -V^# = -v/ P ^, (32) 
dp Jpo dp p 

where p , an integration constant, plays the role of a limiting density. Hence the matter 
Lagrangian can be obtained by a simple integration as 

L m (p)=p[l + U(p)}- ["dp, (33) 

Jpo 

where II (p) = f£ o dp/ p, and we have normalized an arbitrary integration constant to one. 
The quantity p is an integration constant, or a limiting pressure. Therefore the correspond- 
ing energy-momentum tensor of the matter satisfying a barotropic equation of state is given 
by 

T>» = {p [1 + $ (p)} + p (p)}u»u v - p (p) <T , (34) 

respectively, where 

*(P)= f P ^ 2 dp = U(p)-^, (35) 

Jpo p 2 p 

and with all the constant terms included in the definition of p. By introducing the energy 
density of the massive body according to the definition 

e = p[l + $(p)], (36) 
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the energy-momentum tensor of a test fluid can be written in / (R, L m ) gravity in a form 
similar to the standard general relativistic case, 



= [e (p) + V {p)\ - p (p) gT ■ 



(37) 



From a physical point of view $ (p) can be interpreted as the elastic (deformation) poten- 
tial energy of the body, and therefore Eq. ( l34"j) corresponds to the energy-momentum tensor 
of a compressible elastic isotropic system. An energy-momentum tensor of a similar form 



29|]. The matter 



was introduced in the framework of standard general relativity by Fock 
Lagrangian can also be written in the simpler form L m (p) = p$ (p). 

The form of the matter Lagrangian, and the energy-momentum tensor, are strongly 
dependent on the equation of state. For example, if the barotropic equation of state is 
linear, p = (7 — 1) p, 7 = constant, 1 < 7 < 2, then [22J 



L m (p) = p \l + (7 - 1) 



In 



.Po, 



and 



$(p) = ( 7 -l)ln [£. 

\Po 



(38) 



(39) 



respectively. In the case of a polytropic equation of state p = Kp 1+1 ^ n , K,n =constant, we 
find 



L m {p) = P + K 



n 



n + 1 



1 P 



l+l/n 



(40) 



and 



$ (p) = Knp l+l/n = np (p) , (41) 

respectively, where we have taken for simplicity Po = Po = 0. For a fluid satisfying the ideal 
gas equation of state p = kspT/p, where kn is Boltzmann's constant, T is the temperature, 
and p is the mean molecular weight, we obtain 



L m (p) = p \l + 



k B T 



In 



+ Po, 



and 



$(p) 



k B T p 
In — . 

P Po 



(42) 



(43) 



respectively. In the case of a physical system satisfying the ideal gas equation of state, the 
extra-acceleration is given by 



a E 



k B T y 
P 



f Lm {R,L m )\n 



P_ 

Po 



(44) 



and it is proportional to the temperature of the fluid. 

In / (R, L m ) gravity the extra-force does not vanish for any specific choices of the matter 
Lagrangian. In the case of the dust, with p — 0, the extra force is given by 

r = -V, In [f Lm (R, L m )\ (uV - gT) , (45) 

and it is independent on the thermodynamic properties of the system, being completely 
determined by geometry, kinematics and coupling. In the limit of small velocities and weak 
gravitational fields, the extra-acceleration of a dust fluid is given by 

a E = -V[f Lm (R,L m )]. (46) 

If the pressure does not have a thermodynamic or radiative component one can take 
Po = 0. If the pressure is a constant background quantity, independent of the density, so 
that p = po, then L m (p) = p, and the energy-momentum tensor of the matter takes the 
form corresponding to dust, T^ u = pu^u u . 

V. EXPONENTIAL / (R, L m ) GRAVITY 

As a simple toy model for / (R, L m ) gravity we consider that the gravitational field can 
be described by a Lagrangian density of the form 

/ (R, L m ) = A exp (-Lfl +jL m y (47) 

where A > is an arbitrary constant. In the limit (1/2A) R + (1/A) L m 1, we obtain 

f(R 7 L m )^A + ^ + L m + ..., (48) 

that is, we recover the full Hilbert-Einstein gravitational Lagrangian with a cosmological 
constant. 

With this choice of the Lagrangian density the gravitational field equations take the form 

Rfiv — (A — L m ) g^v + T^y — 
1 ^^Vfl + V.VL^^- 



A 



1 

A2 



^V X R + V A /.,„) Qv A /.' + V A /. m ) .y /( „ - 



(49) 
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By introducing the standard Einstein-Hilbert Lagrangian density of the gravitational 
field Lh = R/2 + L m , the gravitational field equations for the exponential / (R, L m ) gravity 
model take the form 



In the case of dust, with L m = p, the mass density of the particles, the equation of motion 
of test particles becomes 



In the case of weak gravitational fields and of small particle velocities, the extra- 
acceleration induced in the exponential model is given by a = — (1/A) [V {R/2) + VL m ]. 
The extra-acceleration induced by the matter-geometry coupling is proportional to the gra- 
dients of the Ricci scalar, and of the matter Lagrangian. In the case of dust, the extra-force 
is proportional to the gradient of the matter density p. 

As one can see from Eq. (|50|) . the gravitational field equations of the exponential 
fi m (R, L m ) gravity model contain a background term, proportional to the metric tensor 
g^p, and which depends on both the constant A, and the physical parameters (density, tem- 
perature, pressure) of the matter. Thus the model introduces an effective, time dependent 
"cosmological constant" . The presence of such a term could explain the late acceleration of 
the Universe, which is suggested by the recent astronomical observations [1]. 

VI. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered a generalized gravity model with a Lagrangian 
density of an arbitrary function of the Lagrange density of the matter, and of the Ricci scalar. 
The proposed action represents the most general extension of the standard Hilbert action for 
the gravitational field, S = f [R/2 + L m ] ^—gd 4 x. Note that the classical Einstein-Hilbert 
action, as well as most of its generalizations, has an additive structure, being constructed 
as the sum of several terms describing independently geometry, matter, and their coupling. 
On the other hand, the / (R, L m ) modified theories of gravity open the possibility of going 




(50) 




(51) 



11 



beyond this algebraic structure - for example in the exponential / (R, L m ) gravity model 
the gravitational Lagrangian density is the product of two independent functions of the 
Ricci scalar and of the matter Lagrange density. The equations of motion corresponding to 
this model show the presence of an extra-force acting on test particles, and the motion is 
generally non-geodesic. The physical implications of such a force have been already analyzed 
in the framework of the generalized gravity model with linear and general coupling between 



matter and geometry, proposed in 



19M2l|. and the possible implications for the dark matter 



problem and for the explanation of the Pioneer anomaly have also been considered. 

On the other hand, the field equations Eqs. are equivalent to the Einstein equations 
of the f(R) model in empty space-time, but differ from them, as well as from standard 
general relativity, in the presence of matter. Therefore the predictions of the present model 
could lead to some major differences, as compared to the predictions of standard general 
relativity, or its extensions ignoring the role of the matter, in several problems of current 
interest, like cosmology, gravitational collapse or the generation of gravitational waves. The 
study of these phenomena may also provide some specific signatures and effects, which could 
distinguish and discriminate between the various theories of modified gravity. 
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